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We consider various generalizations of the Kepler problem to three-dimensional 
sphere S" 3 , a compact space of constant curvature. These generalizations in- 
clude, among other things, addition of a spherical analog of the magnetic 
■ monopole (the Poincare-Appell system) and addition of a more complicated 

. field, which itself is a generalization of the MICZ-system. The mentioned sys- 

tems are integrable — in fact, superintegrable. The latter is due to the vector 
integral, which is analogous to the Laplace-Runge-Lenz vector. We offer a 
classification of the motions and consider a trajectory isomorphism between 
planar and spatial motions. The presented results can be easily extended to 
Lobachevsky space L 3 . 
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^ ; 1 The Kepler problem in 

in \ 

Consider the Kepler problem: a mass point (of unit mass, without losing in general- 
ity) moves in the Newtonian field of a fixed center; the intensity of the gravitational 
interaction 7 is constant. 

In this problem, in addition to the integral of energy 



X 
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Ho = \Y.qt + U (1) 

8=1 

and the vector integral of angular momentum 

M = qxq, (2) 
the equations in three-dimensional Euclidean space M 3 = {g 1; q 2 , q 3 }, 

q\ = U = --, r 2 = q\ + q\ + gf , 7 = const, (3) 

oqi r 

have one more remarkable vector integral, which is due to certain hidden symmetry 
of the Kepler problem. This vector integral is called the Laplace-Runge-Lenz vec- 
tor A = (A l , A 2 , A 3 ). It exists only in the case of Newtonian potential (of all the 
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central potentials) and can be written as follows: 

A = M x q + —q. (4) 

r 

Introducing the momenta p = q, we can rewrite equations (j3J) and integrals (HJ), (0), 
and (j3J) in the canonical form: 

dH . dH 3 

P = ~0 q -> q = ~~dp' p,qew. (5) 

The Poisson brackets for the components of the integrals M and A are 

{M u Mj} = e ljk M k , {Mi, A,} = e ijk A k , {A t , Aj} = -2he ljk M k , (6) 

where h is the constant of energy (JTJ), /i = \p 2 — ^, and e^ k is the Levi-Civita 
symbol. Depending on the value of h, the algebra of integrals © is either so(4) 
(when h < 0) or so(3, 1) (when h > 0). 

Note that, since (714", A) = 0, A is always in the plane of the trajectory. Besides, 
the vector's direction coincides with the direction of the ellipse's major axis, while 
its absolute value is proportional to the eccentricity. 

The algebra of integrals © is an algebra, under which the Kepler problem is 
invariant. Invariance under a global group of transformations (i.e., for example, 
under group 5*0(4) for h < 0) was studied by V. A. Fok G. Gyorgyi [TT] and 
J. Moser jTH]- The latter work contains the most general result, which shows that 
even in the n-dimensional case, the constant energy surface (for h < 0) after suitable 
regularization is topologically equivalent to the bundle of unit vectors tangent to n- 
dimensional sphere S n . 

Note also that the principal dynamical effect of a redundant algebra of inte- 
grals © is the fact that the trajectories of system Q are closed in the configura- 
tional and phase spaces. 



2 The MICZ-system in R 3 . Appell's problem 

Consider one more generalization of the Kepler problem, for which an analog of 
integral (j3J) exists. To this end, in the phase space T*M 3 we specify a noncanonical 
Poisson bracket 

{li, Qj} = 0, Pj} = 5 ip {p i7 p k } = -^e ijk ^ (7) 
and a Hamiltonian 

1 3 2 
ff i = 2& 2 "J + ^ li I 1 — const. (8) 
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Remark. This system (the differential equations of motion) can as well be 
obtained with the standard canonical bracket, but in this case the Hamiltonian would 
contain terms linear in momenta. 

Equations (JTJ), (JHJ) define the MICZ-system (Mcintosh- Cisneros-Zwanziger); it 
describes a particle's motion in the asymptotic field of a self-dual monopole [SJ. It 
was formally studied by Zwanziger (2JJ, Mcintosh and Cisneros ^H] without any 
relevant physical interpretation (see also [7j). 

Consider some special cases of (J7|), ©• The Kepler problem can be obtained 
if we put [i = 0. Putting 7 = and fi = in the Hamiltonian (jHJ, not in the 
bracket (J7|), we have the classical integrable Poincare problem of a particle moving 
in the field of a magnetic monopole. As it was shown by Poincare, the particle's 
trajectories in this case are geodesies of a circular cone. 

P. Appell considered a more general problem of a particle moving in the field 
of a Newtonian center and in the field of a magnetic monopole, assuming that the 
center and the monopole coincide []] . This problem occurs if [i = in (jHJ) (in the 
bracket (JTJ), however, /i / 0). In this case, the trajectory is a conic section in the 
involute of the circular cone, while the integral of areas is preserved during the 
motion. On the cone itself, the trajectories are, generally, not closed. 

There is no analog of the integral A in the Poincare and Appell problems, but 
it exists for the system (JTj), (jHJ) • As for the vector integral of angular momentum M , 
it exists for all the above problems. Indeed j^j, the vector functions 



form the algebra of integrals of (J7J), (JHJ), which is isomorphic to so(4) for H 1 < and 
to so(3, 1) for H 1 > 0. 

Again the trajectories are conic sections, and, since (M , q/r) = — /1, belong 
to the circular cone with cone angle 9 = arccos/x/|A4"| and the axis of symmetry, 
defined by M. 

Various generalizations of the Laplace-Runge-Lenz integrals to dynamical sys- 
tems in the Euclidean space were studied in [To] . 

3 The Kepler problem on three-dimensional sphere S' 
(and on Lobachevsky space L 3 ) 

Consider analogs of the Kepler problem in some simple non-Euclidean spaces of 
constant curvature — three-dimensional sphere S 3 and Lobachevsky space. We will 



M = q x p + n-, 



r 



(9) 




(10) 
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discuss the spherical case in more detail. However, all the results, after appropriate 
revision, can be extended to Lobachevsky space. 

Let the three-dimensional sphere S 3 be embedded into the four- dimensional Eu- 
clidean space R 4 = {g , q±, q 2 , q 3 } and given by equation 

ql + ql + ql + ql = R\ (n) 

where R is the sphere's radius. 

Introduce spherical coordinates on S 3 : 

q = R cos 9, q 1 = Rsin9 cosip, 
q 2 = R sin 9 simp cos ip , q 3 = Rsin9 sincpsinip. 

Consider the motion of a particle in the field of a Newtonian center, placed at 
one of the poles, 9 = 0, of the three-dimensional sphere. It is well known [T71 IT21 
ITU1 IT3*1 ITo^ that the analog of Newtonian potential on sphere is 

U = - 7 cot 9 = - 7 ^, q 2 = q\ + qj + q 2 3 , q = (q lt q 2 , q 3 ). (13) 

Recall that the potential (jlHjl can be obtained either by solving the Laplace- 
Beltrami equation on sphere S 3 (see below, (I29j) ) . This equation is invariant under 
group SO (3) and has a singularity at the pole 9 = 0, or by extending Bertrand's 
theorem to sphere [HI Ej . 

In independent coordinates q = (51,52, 93)) the Lagrangian of the problem in 
question is 

L = ^q 2 + q 2 (q,q) 2 ) 2 -U(q), (14) 

where q is found from namely, q = ±a/ R? — q 2 . After introducing the 

momenta 

9q y/i? 2 - q 2 

the equations of motion can be represented in the canonical Hamiltonian form with 
Hamiltonian 

H=lp 2 -^( P ,q) 2 + V(q). (16) 

These equations have a vector integral of angular momentum 

M = pxq = qxq (17) 

(it exists as well for every "central" potential V that depends only on \ q\ = a/ ' q\ + q 2 + 53) 
and an analog of the Laplace-Runge-Lenz integral 

A = q Q p x M + ^R 2 -^-. (18) 




(19) 



F 1 = (M,A), F 2 




M 2 + (M 2 ) 2 



(20) 



while its symplectic leaf is four- dimensional (i. e., the rank of (JT§j) is four). Here, A = 
1/R 2 is the curvature of the space. For real motions, the Kepler problem gives 



The compactness of the symplectic leaf (}2T|) is defined by the curvature of space, A, 
and the value of the constant of energy: 

1. When A = 0: h < — compact, h > — noncompact. 

2. When A > 0: always compact. 

3. When A < 0: h < and h 2 > •y 2 - - the leaf (|2~T|) is disconnected, one 
component is compact, while the other is noncompact; h > —7 — the leaf is 
connected, but noncompact. 

The trajectories of the Kepler problem on sphere (and pseudosphere) are conic 
sections, the generalization of Kepler's laws to this case was done in E3 IS]- 
In the paper jS], bifurcational analysis of the Kepler problem on S 3 and L 3 was 
performed, and the action- angle variables were introduced (see also 4J). 



4 Generalization of the Poincare and Appell prob- 
lems to S* 3 



First, we obtain the Hamiltonian form of the equations of a particle's motion on 
three-dimensional sphere S 3 under generalized potential forces. Indeed, consider 
the Lagrangian 



where W = W(q) = (W 1 , W 2 , W 3 ) is the vector potential. Introducing the gener- 
alized momenta 



F x = 0, F 2 = 7 2 R\ 



(21) 



L = -{q + q 2 (q,q 



) 2 )-(q,W(q))-U(q) 



(22) 



dL 



q,q 



we obtain the Hamiltonian 

H = i(p + W) 2 - ^(p + W, q) 2 + U(q) (24) 

and the canonical Poisson bracket ({q^pA = S^). Due to a number of consider- 
ations, it is more convenient to study Hamiltonian equations written in terms of 
slightly modified momenta p = p + W, which form the following noncanonical 
Poisson brackets 

{PvP >i 3q 3 % ^' (25) 
{Qi, Pj} = hi, Ui, Qj} = 0, 
where B = rot W. The Hamiltonian (|24j) simplifies, in this case, to: 

H= 1 -p 2 -j- 2 (p,q) 2 + U(q). (26) 

In the case of three-dimensional sphere, an analog of the vector potential of a 
magnetic monopole can be obtained in the following way. 

Electromagnetic field tensor in vacuum satisfies the Maxwell equations 

d a F M + dpFj a + <9 7 F Q/3 = a, p = 0, 1, 2, 3 

d \ (27) 
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where ||<? a /?|| is the metric of space-time, g = det Hs^H- 

For S 3 , the metric of space-time in the spherical coordinates (|12j) is 

dS 2 = c 2 dt 2 - R 2 (d6 2 + sin 2 B[d^ 2 + sin 2 pdif; 2 )) . (28) 

Let i, j, k stand only for the spatial indices, while denotes the spatial portion 
of the metric, taken with the negative sign. We will search for the solution, similar 
to that for a magnetic monopole in flat space, in the form 

F 0l = 0, y/fiF* = e&djj. 

From (|2~7jl we find the equation for the unknown function / 

d k {^gZg tk dJ) = 0, 

coinciding with the Laplace-Beltrami equation. The solution, invariant under group 
SO (3) (i.e. independent of if), <p), satisfies the equation 

* (29) 



sin 
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and looks as follows: 

f = a cot 9, a = const. (30) 

Remark. For Lobachevsky space L 3 , a similar reasoning yields f = acoth^, 
a = const. 

The vector potential of the magnetic monopole W is found from F^ = djWj — 
djW { . In terms of spherical coordinates (9, (p,ip), it reads: 

W e = 0, W v = 0, Wj, = aRcos (p. 

In terms of variables q Q , q, it can be written as 

W=(o, c^^S, -"t^-A^I , (31) 

V W\<il + <iV \q\ <& + <&) 



while for B = rot W, we have 



B = ~q. (32) 



Consider a particle, moving on S* 3 in the field of a Newtonian center and in 
the field of a magnetic monopole, the center and the monopole being placed at the 
pole 6 = 0. This is a spherical analog of the Appell problem. The Hamiltonian of 
the problem is either or (j2HJ) with W(q) and U(q) defined, respectively, by (|3*Tj) 
and (j!3)l . Hamiltonian equations always admit the integral of angular momentum 

M = p x q - a T^7 = q x q - a^-r. (33) 

|g| \q\ 

To simplify the reasoning, we put R = 1 and write the Lagrangian in the spherical 
coordinates (fT2f 

L = - (e 2 + sin 2 B(p 2 + sin 2 9 sin 2 <^ 2 ) + a cos (pip -U{6). (34) 

Suppose that the angular momentum vector (jHHjl is aligned with the axis q 1 in 
the space q x , q 2 , q 3 ; then 

M 2 = Mo = 0, = 0, sin 2 (hji = - = const. (35) 

COS(/9 

The latter relation is a generalization of Kepler's second law. Consider the 

invariant surface in S* 3 given by (^M, j^-J = const. On this surface, choose a point, 

which is 29 away from the particle. Then, the great-circle arc, joining the origin 
of coordinates with the chosen point sweeps equal areas in equal time intervals. 

(Indeed, the time rate of change of the area is ^ = ^ J^ 9 sin o do^j ^ = 2 sin 2 #?/>.) 



Taking into account the integral of energy E = h and (|35j) . we obtain 

e = y/2{h-u e (0)), uM = u(e) + l^4^ (36) 

2 sin 6 

and the explicit expression for the trajectory in terms of quadratures: 

a d6 



s^e x /2(h-u)-^ 



a 2 sin 2 <p 



di/j, (37) 



where h = h cos 2 <£>q, U = U cos 2 ip . 

When 7 = 0, an explicit expression in terms of quadratures for the analog of 
the Poincare problem is obtained from (j37J), and the trajectories are geodesies on 
the invariant cone defined by (M, q/\q\) = const. As it was noted above, in the 
case of the Appell problem's analog the trajectories are conic sections (i. e. ellipses, 
hyperbolas, parabolas) on the plane development of the cone. Generally, these are 
not closed on the cone, but there is a possibility to "adjust" the potential (jTBjl so 
that the trajectories will be always closed in the presence of a monopole. Due to 
this "tuning", the Euclidean MICZ-model can be generalized to the spherical case, 
for which the Laplace-Runge-Lenz integral exists. 

5 Generalized MICZ-model 

Consider a motion in the field of a monopole and in the field with the potential 

U(0) = -j cot9 + -— 7, n = const. (38) 
2 sin 6 

The trajectory in this case is given by (|37|). If 

u = a 2 , (39) 

the trajectory is 



adO 



dip. (40) 



sin 2 # A /2/i + 27 cot# 



The trajectory (jlUj) is closed, and the gnomonic projection gives us the conic 
section 

tan0 = - V — t pr (41) 

1 + ecos(^ — ipo) 

with the following focal parameter and eccentricity: 



a 2 2a 2 /, a 2 

P= 9.. , e = \ l + .,o 9.. k 



7 cos 2 ip ' y 7 2 cos 2 (p \ 2 cos 2 <p J 
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The expression in terms of quadratures for 9 (I36|) is 



9 2 = 2h + 2 7 cot 9 



a 



2 



cos 2 v? • sin 2 9 
c 2 



= 2/1 + 27 cot 9 ^- = f(9, c, h), 

sin 9 

where c = a 2 /cos 2 y? . Let us plot the bifurcational diagram of the problem's solu- 
tions on the parameter plane (c 2 , h). To this end, recall that at critical points (c + , h t ) 
on a bifurcation curve, f(9 , c^, aj = }'q{9 0) c^, aj =0. As a result, we have two 
curves (Fig. 1): 

I. 2h = c 2 - \; II. c 2 = 0. 

c 2 

Besides, since c = a 2 / cos 2 <p , the inequality c 2 > a 4 also holds. Thus, in the plane 
defined by the constants h, c 2 (see Fig. 1), the domain of allowable values h, c 2 lies 
above the line c 2 = a 4 and below the hyperbola defined by I. For a point of this 
plane above h = ^c 2 , the particle moves only in the upper half-plane of the sphere. 
Otherwise, the particle can also move into the other half-plane. 

It is easy to formulate an analog of Kepler's third law, coinciding with the tra- 
ditional law for a curved space |13J. Indeed, since sin 2 9ip = c, we have 

sin 2 9dib 
dt = -. 



Therefore, 



c Jo c J_ n p 2 + (I + ecosip) 



7T 




(43) 



This dependence of the orbital period on energy can be easily transformed into the 
dependence on (angular) length of the semi-major axis 



T= — \/- tana + + tan 2 a(l + tan 2 a)" 1/2 , (44) 



where tana = — ?. 

h 



As in the Kepler problem (on R and S ), closedness of the trajectories is closely 
connected with some hidden symmetry of the problem, i.e., with existence of a 
vector integral of the Laplace-Runge-Lenz type. 
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For the system (|38|) . ()39|1 this vector can be written as 



A = q p x M + ^R 2 -^. (45) 

The Poisson brackets for the components of A and the components of the integral 
of angular momentum ()33|) are: 

{M i5 Mj} = -e ijk M k , {M t , R 3 ] = -e ijk R k , 



{i^ ^} = 2e yfc [R 2 h — M + M fc . 

As before, we can specify the conditions (in terms of the curvature of the space 
and the value of the integral of energy), under which the symplectic leaf ()23|) is 
compact. 



6 Trajectory isomorphism for central potential sys- 
tems on S 2 and M 2 

For U — U(r), the equations © define a central potential system on R 3 . If U — U (8) 
in the Lagrangian (|14jl. then we have a central potential system on S* 3 . These 
systems, respectively, have flat (M 2 ) and spherical (S 2 ) invariant manifolds. These 
two-dimensional systems can be shown to be related, using the central (gnomonic) 
projection (from the center of the sphere tangent to the plane at the attracting 
center) and some suitable change of time. 

Following Serret and Appell 12] , consider a system in M. 2 with the following 
equations of motion (in polar coordinates): 

dt\dp) ' dt\dv) ' [ ' 

where T p is the kinetic energy of a point in the plane, 

T p =^(p 2 + pV), (48) 

while R, $ stand for certain generalized forces (generally, non-potential). 

Let us perform the transformation of coordinates (the gnomonic projection from 
Fig. 2), forces and time: 



p = tan 9, ip = ip, dt = cos 2 Qdr 

R = cos 2 ee, <s> = cos 2 ey. 



2 n r\ /f> „„„2 n iTr 

(49) 
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This results in a system on S . 



d fdT s \ d fdT t 



Q > 3Z Ul7 = ( 5 °) 



dr \d9' J ' dr \dil)' 
where 0' = ^p-, ■?/>' = while T s is the kinetic energy of a point on the sphere, 

T s = \ {ff + sin 2 ei,' 2 ) • (51) 

It is easy to see that 

Statement. There exists a trajectory isomorphism between the Lagrangian sys- 
tem in R 2 ; with central potential 

L= l -(p 2 + pH 2 ) + U(p), 

and the Lagrangian system on S 2 , with central potential of the form 

L = l(fl 2 + sin 2 dip 2 ) + U {tan 6). 

To prove that, it is sufficient to put in (J47|) - (|51j) 

$ = ^ = o 



dU dU dU dp R 

R = — ^ — i <y = — 



dp ' 89 dp 86 cos 2 9 

These transformations easily bring the plane Kepler problem to its analog on 
sphere. 

Note that under the transformation ()49|) a potential field of forces can be trans- 
formed into a non-potential field, and vice versa. 

If we consider the inverse transformation (J49|) as a transformation from sphere 
to plane, we have to adopt negative values of p. In this case, p is negative when 
the trajectory on the sphere crosses its equator, or, in the plane (p, ip), when the 
trajectory jumps from +oo to — oo. If, instead of p = tan. 9, we consider £ = cot#, 
then the trajectory in the plane (£, ip) is continuous. 

It is also easy to show that the described isomorphism can be extended to the 
generalized potential systems discussed in Sections IH El Note also that the trans- 
formation (|49p. applied to the Kepler problem on S 2 , can be used to generalize the 
Bohlin (Levi-Civita) regularization. It can be shown that on a fixed energy level, 
the Kepler problem is reduced to the harmonic oscillator problem. 
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Another interesting property of the Kepler system in R 2 (due to Hamilton) is 
that the velocity hodograph for a moving point is a circle with a displaced center. 
A similar vector can be specified for the Kepler problem on S 2 : 

x 

n= l + x 2 /R 2 ' 

where x = (p cos (p, p sin <p) is the radius- vector of a point under the gnomonic 
projection ([49)1 . The tz hodograph can be easily shown to be a circle with a displaced 
center. 
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Figure 1: 




Figure 2: The gnomonic projection 
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